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Abstract. We show how the magneto-phonon resonance, particularly pronounced in 
sp 2 carbon allotropes, can be used as a tool to probe the band structure of multilayer 
graphene specimens. Even when electronic excitations cannot be directly observed, 
their coupling to the E 2g phonon leads to pronounced oscillations of the phonon 
. feature observed through Raman scattering experiments with multiple periods and 

I amplitudes detemined by the electronic excitation spectrum. Such experiment and 

analysis have been performed up to 28T on an exfoliated 4-layer graphene specimen 
deposited on Si02, and the observed oscillations correspond to the specific AB stacked 
4-layer graphene electronic excitation spectrum. 
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1. Introduction 

Since the theoretical prediction of magneto-phonon resonance in graphene and bilayer 
graphene [H [2j [3], the resonant interaction between inter-Landau level electronic 
excitations and T-point optical phonon in graphene, this effect has been widely explored 
experimentally in neutral graphene-like systems such as multi layer epitaxial graphene 
on C-face SiC [4] , graphene like inclusions on the surface of bulk graphite [5l [6] , and 
in bulk graphite pD, E]. Magneto-phonon resonance in these systems are pronounced 
because (i) they are gapless and low energy direct (Ak = 0) electronic excitations exist, 
and (ii) because of the Kohn anomaly [9] [H)] at the T point of the phonon Brillouin 
zone which makes phonon energies ultra sensitive to modifications of the electronic 
excitation spectrum. Such modifications can be achieved by tuning the Fermi energy or 
by applying a magnetic field perpendicular to the plane of the quasi-2D crystal. 

Performed in systems with a well-known electronic excitation spectrum, magneto- 
phonon resonance experiments provide precise information concerning the electron- 
phonon interaction through the amplitude of the observed oscillations of the energy and 
of the line width of the phonon feature observed in Raman scattering experiments, as a 
function of the applied magnetic field. On the other hand, since electronic excitations 
which couple to the E 2g optical phonons are well identified by symmetry arguments 
(optical-like excitations characterized by a change of angular momentum A|n| = ±1 pQ, 
where n is the Landau level index), this effect can be used to determine the electronic 
excitation spectrum. The evolution of optical-like excitations, calculated in the frame 
of Ref. [11], as a function of the magnetic field is presented in figured] for a) graphene, 
b) bilayer graphene, c) tri-layer graphene and d) for four-layer graphene. Each of these 
distinct electronic excitation spectra are tuned in resonance with the phonon energy at 
different characteristic values of the magnetic field. In analogy to magneto-transport 
experiments which probe the density of states at the Fermi energy, magneto-phonon 
resonance can be used as a probe of A|n| = ±1 electronic excitation at the E 2g phonon 
energy, and Landau level spectroscopy can be achieved by modelling the observed shifts 
of the phonon feature. 

In this paper, we report on low-temperature magneto-Raman scattering 
experiments revealing the magneto-phonon resonance of an exfoliated multi-layer 
graphene specimen deposited on a Si0 2 /Si substrate. The investigated sample is 
characterized by a low "intrinsic" doping and oscillations of both the phonon energy 
and line width can be observed for magnetic fields B > 7 T. These oscillations have a 
noticeably more complicated pattern than for monolayer graphene and for bulk graphite. 
From the analysis of this pattern we deduce the number of layers in the sample, which 
turns out to be equal to four. Observation of a magneto-phonon resonance can hence 
be used to trace electronic excitations at the phonon energy and determine the band 
structure of the probed electronic system. 
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Figure 1. Evolution of inter Landau level excitations with A|n| — ±1 calculated in 
the frame of the model developed in Ref. |TT| , in the case of (a) a monolayer graphene, 
(b) a bilayer graphene, (c) a trilayer graphene and (d) four-layer graphene. The values 
of k z used in these calculations are indicated in the figures. The phonon energy is 
indicated by the solid horizontal line at 1585 cm" 1 . Dashed horizontal lines indicate 
resonant magnetic fields corresponding to the different excitation spectra. 

2. Experimental details 

The multilayer graphene specimen is presented in the inset of figure |2k. It has been 
exfoliated from bulk graphite and deposited on a Si/Si0 2 substrate with 300 nm of 
Si0 2 . A single metallic electrode was then deposited on the flake (not used in this 
experiment). To locate the sample, we spatially map the Raman scattering signal of 
the silicon substrate, which shows a strong attenuation when the laser spot is focused 
on the metallic layer. A typical room temperature Raman scattering spectrum of the 
multilayer graphene flake together with an optical image are presented in figure EK The 
2D band has a multi-components broad feature. We find the standard analysis of 2D 
band features in multi-layers specimens, that is, the deconvolution of the complex 2D 
band line shape using 2N Lorentzian functions (where N > 1 is the number of layers) , 
relevant for mono-layers, for bi-layers and for tri-layers [T2l [TBI H%] . not precise enough 
for N > 3. As a consequence, 4-layer graphene flakes can hardly be distinguished 
from the 5-layer and additional information is needed to fully characterized multilayer 
graphene systems. 




Figure 2. a) Raman scattering spectrum of the four-layer graphene specimen 
measured at room temperature with the A = 514.5 nm line of an argon laser. Inset: 
Microscope image of the 4-layer graphene flake, b) Schematics of the micro-Raman 
scattering set-up operating at liquid helium temperature and in high magnetic fields. 



The low-temperature magneto-Raman scattering response of our multilayer 
graphene has been measured with a home-made micro- Raman scattering setup presented 
in figure [2b. Laser excitation at A = 488 nm produced by an Ar+ laser is injected into 
a 5 /im core mono-mode optical fiber which is fixed on a miniaturized optical bench and 
placed at the center of the magnet. This experimental setup can host various optical 
filters which are used to filter the excitation laser light at the output of the mono 
mode fiber (laser line filter), to selectively reflect the excitation laser and transmit the 
scattered signal (dichroic mirror), to reject the scattered laser before the collection fiber 
(notch filter), and to impose a circular polarization (linear polarizers and quarter wave 
plate). The excitation spot on the surface of the sample has a diameter of ~ 1 /im and 
the optical power is of 5 mW. The sample is mounted on piezzo stages which allow 
to move the sample under the laser spot with a sub-micrometer resolution. Polarized 
Raman scattering spectra were measured in nearly back-scattering Faraday geometry. 
The collected light was dispersed with a single grating spectrometer (spectral resolution 
AA = 0.07 nm) equipped with nitrogen cooled CCD detector and band pass filters were 
used to reject the stray light. 

3. Magneto-Raman scattering of exfoliated multilayer graphene specimen 

In order to obtain further information concerning our multilayer graphene sample, we 
have performed magneto- Raman scattering experiments in the G-band range of energy. 
The G-band feature of the Raman scattering spectrum is observable in the crossed 
circular polarization configuration [6] which is used in the present experiment. Even 
though electronic excitations cannot be, in the present case, directly imaged on the 
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Figure 3. False color map of the scattered intensity in the energy range of the E^g 
phonon as a function of the magnetic field in the o~ /a + (a) and er + /er~ (b) polarization 
configurations. 



Raman scattering spectra, their effect when tuned in resonance to the E 2fl optical phonon 
by applying a magnetic field can be clearly observed. 

The magneto-phonon resonance of our multilayer graphene sample is presented 
in figure |3] in the form of false color map of the scattered intensity as a function of 
the magnetic field, in the two crossed circular polarization configurations <T ± /cr T . the 
characteristic individual spectra in a~/a + configuration are presented in figure UK In 
both configurations, we observe, for B > 7 T, oscillations of the energy of the Ei g 
phonon and of its line width. This low value of the magnetic field for the onset of 
magneto-oscillations is an indication that the carrier density in this sample is lower 
than half of the Ei g phonon energy. The resonant values of the magnetic field are 
different in the two polarization configurations, B res = 7, 7.85, 10, 11.9, 16.75, 26.6 T 
in the cr + / ' o~ configuration and B res = 6.7, 7.3, 9.7, 11.4, 15.7, 25.6 T in the a~/a + 
configuration. Moreover, the electron-phonon coupling appears to be more efficient at 
particular values of the magnetic field (indicated above by the bold typeface) where 
the amplitude of the observed oscillation is enhanced with respect to other resonant 
magnetic fields. This nontrivial behavior is different from the one observed for epitaxial 
graphene, for graphene-like domains on the surface of bulk graphite or for K-point 




Figure 4. a) Raman scattering spectra measured in the a~/a + polarization 
configuration at different values of the magnetic field. The dashed vertical line indicates 
the B=0 energy of the Ei g phonon. b) Intensity of the scattered light at the phonon 
energy at B=0 as a function 1/B. c) Amplitude of the Fourier transform of b) showing 
two characteristic frequencies. 



carriers in bulk graphite, where the amplitude of the oscillations was found to increase 
with the magnetic field as the effective electron-phonon coupling strength depends on 
the degeneracy of the Landau levels which increases linearly with B. 

To determine the number of Landau level fan charts contributing to the observed 
magneto-phonon resonance, we extract from the spectra presented in figure [3^ and in 
figure the evolution of the amplitude of the scattered light at the phonon energy at 
B = as a function of the inverse magnetic field. This quantity shows oscillations which 
are representative of the phonon oscillations: each time the phonon feature splits or is 
shifted due to the resonant interaction with an electronic excitation, the amplitude of 
the scattered light at the phonon energy at B = decreases. This evolution is presented 
in figure Hb. It is then possible to numerically calculate the Fourier transform of this 
signal, presented in figure Hfc, which shows two distinct frequencies at about 21 and 
49 T, corresponding to two periods A(l/B) of about 0.05 T _1 and 0.02 T _1 . 

Their appearance can be understood qualitatively using the results of Ref. [H] 
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for the electronic spectrum of a iV-layer graphene. Namely, the spectrum is similar to 
that of a bulk graphite, but instead of the continuous quantum number k z , —tt/2 < 
k z < tt/2 (the wave vector, perpendicular to the layers, in the units of the inverse 
interlayer spacing), only N/2 (for N even) or (N + l)/2 (for iV odd) discrete values, 
k z = n/(N+l), . . . , (N/2)ir/(N+l) or k z = n/(N+l), 2n/(N+l), . . . , tt/2, respectively, 
are allowed. Now, instead of a continuum of inter-Landau-level excitations for all k z in 
graphite, which couple to the phonon [7], in the 4-layer graphene one has two series of 
inter-Landau-level excitations, corresponding to the two allowed values of k z . Thus, the 
presence of two periods indicated that our sample has either three or four layers. 

In the two-band approximation for the effective bilayer model of Ref. [11] . the values 
of the magnetic field B n at which the magnetophonon resonance occurs, are determined 
by 



n(n - 1) + Jn(n + 1)) = u ph , (1) 

v / 71 cos k z 

where Vp is the Fermi velocity in the monolayer, and 71 is the nearest-neighbor inter- 



layer matrix element. Approximating ^jn{n ± 1) w n ± 1/2, we obtain 
1 2vle 



, n. (2) 

B n W p h7i cos k z 

The coefficient in front of n thus gives the period A(l/B) for each value of k z . For the 
monolayer- type spectrum at k z = tt/2, the same consideration (again, approximating 
n(n ± 1) ~ n ± 1/2) gives 

1 8v 2 F {e/c) 



Bn U^h 



(n + 1/2). (3) 



For cjph = 196 meV and vp = 1.01 x 10 6 m/s this gives an inverse period of 7.1 T, which 
is too small compared to the smallest observed 21 T. Thus, we exclude the possibility 
that our sample could be a trilayer graphene. For four layers, we have k z = 7r/5,27r/5, 
which, for 71 = 380 meV, gives two frequencies of 16 and 41 T. This is an indication 
that the sample has four layers. 

To give more solid arguments, we perform a quantitative modelling of the spectra. 
The Raman scattering spectra presented in figure H] can be described by Lorentzian 
functions to extract the energy corresponding to the maximum of scattered light and the 
full width at half maximum (FWHM) of the phonon feature. These results are presented 
by the black dots in figure [5] (a-d) for both polarization configurations. Again, the 
observed behavior is not a simple series of oscillations; rather, two series of oscillations 
can be distinguished. Qualitatively, this is explained by the same arguments as above 
(two allowed values of k z ). The theoretical curves are shown in figure. [5] as solid curves 
on top of the experimental symbols. 
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Figure 5. Evolution of the maximum of scattered light and of the FWHM of the G 
band as a function of the magnetic field in the o~ /a + (a and b) and in the a + /cr~ 
polarization configurations. Red solid lines are the result of the calculation detailed in 
the text. 



4. Theoretical aspects of the magneto-phonon resonance in 4-layer 
graphene 

In contrast to bulk graphite, where due to the translational invariance in the z direction, 
the momentum conservation holds, which implies that only phonons with the wave 
vector q z = are probed (assuming the photon momentum to be very small), and 
that the electronic transitions involved are those from states with some k z into states 
with k' z = k z — q z = k z , the 4-layer graphene lacks the translational invariance in the 
z direction, which introduces some additional difficulties in the theoretical description 
of the magneto-phonon resonance, as compared to the bulk graphite. 

(i) There is a priori no reason to focus on one particular phonon mode. In total, 
there are eight degenerate optical phonon modes with zero in-plane wave vector (two 
for each layer). The circular basis for the in-plane polarizations is still a good one, 
thus for each circular polarization there are four phonon modes. We assume that their 
degeneracy is lifted by the electron-phonon interaction, neglecting the elastic coupling 
between the layers. Thus, for each of the two circular polarizations, the normal modes 
should be found as eigenvectors of the electronic polarization operator, which is a 4 x 4 
matrix. 

(ii) The electronic polarization operator not only has a contribution from the k z - 
diagonal electronic transitions, but also from transitions between states with different k z , 
since no momentum conservation holds. This gives not just two series of transitions, 
but four, corresponding to all possible combinations of {n/5,27r/5} — > {iy/5,2ti/5}. 
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Still, it can be naturally expected (and is confirmed by the actual calculation) that the 
overlap between states with different k z is much smaller than that between states with 
the same k z , this is why two series of oscillations instead of four are seen in figure HI 

(iii) Generally speaking, all four phonon modes should contribute to the Raman 
spectrum. The intensities of the four components are determined by the projections 
of the corresponding mode eigenvectors on the Raman matrix element (photon- 
photon-phonon vertex), which is also a four-component vector for each given circular 
polarization configuration. In bulk graphite, due to the momentum conservation, the 
Raman vertex would be simply proportional to a 5-function 5(q z ), but in the 4-layer 
sample it has a non-trivial structure. 

The electronic band structure of multilayer graphene is modelled similarly to 
Koshino and Ando Namely, in the nearest-layer tight-binding approximation, the 
band structure of any Bernal-stacked multi-layer graphene can be deduced from the 
one of the bulk graphite by selecting specific values of k z corresponding to the possible 
standing waves perpendicular to the basal plane of the multilayer graphene specimen. 
The monolayer graphene then corresponds to a single value of k z = tt/2, the bilayer to 
k z = 7r/3, the trilayer corresponds to two values of k z = 7r/4, tt/2, etc. At each value of 
k z , the Hamiltonian has the form of an effective graphene bilayer: 



r 2 

v F p + 




v F P- 

r 5 



-a^v F p. 
v F p + 





v F p- 

r 2 



(4) 



where 



C = cosfc. 



v F 



3 70a 

2 h ' 



r 1 = 2 7l c, r 2 = 2 72 c 5 
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7o 



c, r 5 = 2 75 C 2 +A, (5) 



a = 1.42 A is the distance between the neighboring carbon atoms in the 
same layer, p± = —ih(d x ± id y ) are the in-plane momenta counted from the corners 
of the hexagonal first Brillouin zone, and 70 = 3.1 eV, 71 = 0.38 eV, 7 2 = —20 meV, 
74 = 44 meV, 75 = 94 meV, A = —8 meV are the parameters of the tight-binding 
Slonczewski-Weiss-McClure (SWM) model [TH [16], whose values we take the same as 
in Ref. [7j. As in Ref. [7], we have neglected the 73 term, responsible for the trigonal 
warping, in order to simplify the diagonalization of the electronic Hamiltonian with an 
external magnetic field. The 7 2 ,7s terms in the SWM model correspond to the second- 
nearest-layer coupling, and their non zero values modify the electronic wave functions 
on the layers with respect to the simple recipe of Ref. [IT] . Here, we consider the 
wave functions of Ref. pj] but we still include the 7 2 and 75 terms in Eq. (jlj), since 
the experimental transition energies are better reproduced when also considering these 
terms. Thus, even though Eq. (jlj), strictly speaking, does not correspond to the SWM 
model, we use it as an effective electronic Hamiltonian for the 4-layer graphene. Given 
Eq. (HI), the calculation of the polarization operator is analogous to the one presented 
in Ref. [7j with summation over k z , k' z instead of integration over k z . Following points 
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(i) and (ii) discussed at the beginning of this section, we also introduce a layer index 
for the electron-phonon coupling. The solid curves presented in figure. |5] were obtained 
using the values A = 4 x 10~ 3 , 5 = 100 meV for the dimensionless electron-phonon 
coupling constant and the electronic broadening respectively, but the values A = 3x 10~ 3 , 
5 = 80 meV or A = 5 x 10~ 3 , 5 = 120 meV differ very little. (We remind that the main 
purpose of the present modelling is not to extract precise values of parameters but to 
demonstrate that the sample is quadri-layer graphene.) 

The third difficulty, mentioned at the beginning of this section, can be bypassed 
since we are interested only in the positions and widths of the Raman peaks (determined 
by the complex eigenvalues of the electronic polarization operator), and not in their 
intensities. Moreover, it turns out that the four eigenvalues behave quite differently 
upon varying the magnetic field. Only one eigenvalue exhibits pronounced magneto- 
oscillations, while other three remain close to the bare phonon frequency. This behavior 
can be understood by noting that near each resonance a, the polarization operator can 
be approximated as IIy(u;) ~ vfiy'*)*/ (u — u a + iT a ), where u a and T a are the frequency 
and the broadening of the corresponding transition, and vf is the corresponding matrix 
element of the electron-phonon coupling (the indices i, j = 1,2,3,4 label the layers where 
the phonon displacements occur). Then, near the resonance, one phonon eigenvector is 
along vf, and other three are orthogonal to it, so that they are effectively decoupled from 
the resonance. The theoretical curves in figure [5] correspond to the eigenvalue exhibiting 
the non-trivial oscillating behavior. Nevertheless, the contribution of other modes, not 
subject to oscillations, can also be distinguished in the spectrum of Fig. |3k) around 
the strongest resonance at B = 25 T. Finally, we note that the possibility to avoid 
the calculation of the Raman matrix element in the 4-layer structure is quite fortunate, 
since its calculation requires a significant effort even in the monolayer graphene |17j . 

5. Conclusions 

We have shown that magneto-phonon resonance in multi-layer graphene specimens can 
be used to determine the specific band structure of these systems. We demonstrate 
this possibility in this work with a 4-layer graphene flake deposited on SiC>2. Magneto- 
phonon resonance in this sample result in a series of oscillations of the phonon energy 
which can be directly related to resonances between two distinct Landau Levels fan 
charts, characteristic of the 4-layer graphene system. 
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